Source breakup dynamics in Au+Au Collisions at y / s/v r /v=200 GeV via 
three-dimensional two-pion source imaging 
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A three-dimensional (3D) correlation function obtained from mid- rapidity, low px, pion pairs in 
central Au+Au collisions at y / ijvjv=200 GeV is studied. The extracted model-independent source 
function indicates a long range tail in the directions of the pion pair transverse momentum (out) 
and the beam (long). Model comparisons to these distensions indicate a proper breakup time 
to ~ 9 fm/c and a mean proper emission duration At ~ 2 fm/c, leading to sizable emission time 
differences ((|AtLCM|) ~ 12 fm/c) partly due to resonance decays. They also suggest an outside-in 
"burning" of the emission source reminiscent of many hydrodynamical models. 
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Collisions between heavy nuclei at ultra-relativistic en- 
ergies produce transient systems with energy densities 
much greater than that required to decompose bulk nu- 
clear matter into quarks and gluons [lj]. Such systems 
were predicted to have long lifetimes if a first order phase 
transition occurred during their formation or decay 0]. 

A number of interferometry studies Q have been made 
to search for signals of long time delays in emissions from 
actual reaction sources [4|. For a Gaussian source func- 
tion, assumed in the traditional Hanbury Brown Twiss 
(HBT) methodology, this would be signaled by an in- 
crease in the width R of the emission source function 
in the out direction of the Bertsch-Pratt coordinate sys- 
tem i.e -Rout/^sidc >> 1- No such result has been 
found by these HBT studies, and the reported Gaus- 
sian source functions are spheroidal with R out « -R s ide in 
the longitudinally co- moving system (LCMS) How- 
ever, a recent study with a ID source imaging tech- 
nique 0, @, 0j nas observed a long non-Gaussian tail 
in the radial source function and attributed it to possible 
lifetime effects [H, [t| . This suggests that further study of 
the source image may give new insights into the reaction 
dynamics leading to source breakup. 

Here, we extract and perform a detailed study of the 
3D two-pion source function using the technique pro- 
posed by Danielewicz and Pratt [lCj, [ll|. Namely, the 
3D correlation function is first decomposed into a basis 
of Cartesian surface-spherical harmonics to extract the 
coefficients, also called moments, of the expansion. In 
turn, they are then imaged or fitted with a trial function 
to extract the 3D source function, which is then used 
to probe the emission dynamics of the two-pion source 
produced. 

Au+Au data (at v /sjvjv=200 GeV) was recorded dur- 
ing 2004 with the PHENIX detector [12] at the Rela- 
tivistic Heavy Ion Collider (RHIC). The collision vertex 
z (along the beam axis) was constrained to \z\ < 30 cm 
of the nominal crossing point. Charged pions were de- 
tected in the east and west central arms of PHENIX, 
each of which subtends 90° in azimuth <f), and ±0.35 units 
of pseudo-rapidity r\. Tracking and momentum recon- 
struction were accomplished with the drift chamber and 
two layers of multi-wire proportional chambers with pad 
readout (PCI and PC3). Particle momenta were mea- 
sured with a resolution cr p /p = 0.7% 0l.O%p (GeV/c). 

Pion identification was achieved for pt < 2.0 GeV/c 
and pt < 1 GeV/c in the time of flight and 
electromagnetic calorimeter, respectively. For this 
analysis, mid-rapidity pion pairs were selected with 
0.2 < p T < 0.36 GeV/c from semi-central (0 - 20%) 
Au+Au collisions at v / J/vtv=200 GeV, where pr is half 
the pion pair transverse momentum. Track merging and 
splitting effects were removed by appropriate cuts on 
both the real and mixed pair distributions System- 
atic variations of these cuts were explored to obtain sys- 
tematic error estimates; on average, they are well within 



the statistical uncertainty. Hence, the pair cuts do not 
introduce any significant bias in the correlation function. 

The 3D correlation function C(q) = -/Vf g d(q)/Abkg(q) 
was constructed as a ratio of 3D relative momentum dis- 
tribution for 7r + 7r + and 7r~7r~ pairs in the same event 
-^fgd(q) to that from mixed events -/Vbk g (q)- C(q) is 
normalised to unity for 50 < |q| < 100 MeV/c. Here, 
q = ( pi ~ P2 ) where pi and P2 arc the momentum 4- 
vectors in the pair center of mass system (PCMS). The 
Lorentz transformation of q from the laboratory frame to 
the PCMS is made by a transformation to the pair LCMS 
along the beam direction followed by a transformation to 
the PCMS along the pair transverse momentum [l3j |. 

To obtain the moments, the 3D correlation function 
C(q), is expanded in a Cartesian surface-spherical har- 



monic basis 



C(q) - 1 = i?(q) = ]T E <...«, (8) <...a,(fiq)> 

/ ai...ai 

(1) 

where I = 0,1,2,..., a t = x,y or z, A l ai ai (f2 q ) are 
Cartesian harmonic basis elements; (J7 q is the solid an- 
gle in q space); R l a a , (q) are Cartesian correlation mo- 
ments given by Eq. ©; and q is the modulus of q. 



(2l + l)\\ fcM^, 
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j^^L.-a^m- (2) 



Here, the coordinate axes are oriented so that z (long) 
is parallel to the beam direction, x (out) points in the 
direction of the total momentum of the pair in the LCMS. 

Correlation moments, for each order I, can be calcu- 
lated from the measured 3D correlation function using 
Eq. ((2]). Eq. (fl]) is truncated at / = 6 and expressed in 
terms of independent even moments only. As expected 
from symmetry considerations, odd moments were found 
to be consistent with zero within statistical uncertai nty ; 
higher order moments were found to be negligible [14j . 
Up to order 6, there are 10 independent moments: 



Rx2 1 Ry2> Rx4i Ryii R'x2y2i ^a6) Ry6' ^x4y2 an d ^x 

1^2 is shorthand for R xx etc (the dependent mo 



x x2y4 



where R x ^ 

ments are obtained from the independent ones 1^, lit). 
These independent moments were extracted as a function 
of q, by fitting the truncated series to the measured 3D 
correlation function with the moments as the parameters 
of the fit. 

Figure Q] shows the correlation moments R l a a , up to 
order 1 = 6. In panel (a), R°(q) is shown along with 
R(q) = C(q) — 1; both represent angle-averaged corre- 
lation functions, but they are obtained by very different 
methods. For R°(q) one uses Eq. ([2]) and the 3D corre- 
lation function. For R(q), one simply evaluates the ID 
correlation function directly as in Ref. @. 

The very good agreement between R°(q) and R(q) un- 
derlines the absence of any significant angular acceptance 
issues and attests to the reliability of the moment extrac- 
tion technique used. Figures QJb)-(j) show results for the 
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FIG. 1: Experimental correlation moments R l (q) for I =0, 2, 
4, 6. Panel (a) also shows a comparison between R°(q) and 
R(q). Systematic errors are less than the statistical errors. 
The solid lines indicate the fit with the Hump function Eq. (0 . 



moments 1= 2, 4 and 6. Contributions decrease with in- 
creasing I in each direction and are relatively small for 
I =6. This justifies the truncation of the series Eq. ^ 
at I =6. The 3D source function S(r) is obtained from 
these moments via imaging or fitting as discussed below. 

Analogous to Eq. (fTJ), S(r) can be expanded in a Carte- 
sian Surface-spherical harmonic basis (Eq. ([3])) 

5 W=E E S l ai ... ai (r)A l ai ^ ai (Q r ). (3) 

l ax...ai 

Substitution of the series for i?(q) and S(r) into the 3D 
Koonin-Pratt equation 



gives a set of ID relations (Eq. ((5])) [ToL [ll | 

R l ai ...aM) = ^ I dr^K^Si^ir), (5) 



which connects the correlation moments R l ai (q) and 
source moments S a Qi ( r )- ^(r) gi ves the probability 
of emitting a pair of particles with a separation vector 
r in the PCMS. The 3D Kernel, if(q, r), incorporates 
Coulomb interaction and Bose-Einstein symmetrization. 
Strong interaction is assumed to be negligible forpions. 
The ID imaging code of Brown and Danielewicz [B[, @, 
was used to numerically invert each correlation moment 
B? a a (q) to extract the corresponding source moment 
S l a „ a (r); the latter were then combined as in Eq. (J3]) 
to obtain the source function. 

The 3D source function can also be extracted by di- 
rectly fitting the 3D correlation function with an assumed 
functional form for S(r). This corresponds to a simulta- 
neous fit of the ten independent moments. A 4-parameter 
3D Gaussian (ellipsoid) fit, using MINUIT minimization, 
gives a poor result (x 2 /ndf=3.7). The solid line in Fig. Q] 
shows the result of a fit to the independent moments with 
an empirical Hump function given by 



S H (r x ,r y ,r z ) = Aexp[-/ S (-^- + + -^-) 



-M 



Ar 2 



4r£ 



4r 2 

ys 



4r 2 

y 2 z 2 



)], (6) 



C(q) - 1 = i?(q) = / dvKici, r)5(r), (4) 



where A, ro, r xs , r ys , r zs , r x \, r y \, r z \ are fit parameters and 
j s = !/[! + {r/r ) 2 ], f\ = 1 - / s . This 8-parameter Hump 
function achieves a better fit to the data (% 2 /ndf=1.4). 
Smearing the track momenta by the measured resolution 
has a negligible effect on the data points and fits. 

Figure [U[a)-(c) shows a comparison of profiles of the 
two-pion 3D source function in the cc, y and z direc- 
tions (S(ra:), S(r y ) and S(r z )) obtained via fitting (line) 
and source imaging (squares). Source image extraction 
makes no assumption for the shape of the 3D source func- 
tion; whereas, moment fitting explicitly assumes a shape. 
Therefore, the good agreement from the two extraction 
methods confirms the sufficiency of the Hump function 
but not its uniqueness. 

The function S(r x ) is characterized by a long tail, 
which is resolved up to ~60 fm, in contrast to S(r y ) and 
S(r z ) which range up to ^25 fm. This difference is also 
reflected in the respective correlation profiles (Fig. [2^d)- 
(f)) obtained by summation of the data (circle), fit (line) 
and image (square) moments up to order / = 6 (Coulomb 
effects are not removed). The broader S(r x ) is associated 
with the narrower C(q x ) (Fig.f^a) and (d)), as expected. 

The extended tail lies along the total momentum of 
the pair in the LCMS. Thus, the relative emission times 
between the pions (including those from resonances), as 
well as the source geometry, will contribute directly to 
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FIG. 2: Source function profiles S^z), S(r y ) and S(r z ) (left 
panels) and their associated correlation profiles C(q x ), C(q y ) 
and C(q z ) (right panels) in the PCMS. Symbols are as indi- 
cated. The bands indicate statistical and systematic errors. 
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FIG. 3: Source function comparison between Therminator 
calculation and image for (a) S(r x ), (b) S(r y ), (c) S(r z ) in 
PCMS. Panel (d) compares AtLCM from Therminator events 
with various assumptions for At and resonance emission. 



S(r x ). The lifetime of the source contributes to the range 
of S(r z ), and S(r y ) reflects its mean transverse geomet- 
ric size. The difference between S(r x ) and S(r y ) is thus 
driven by the combination of the emission time difference, 
the freeze-out dynamics and the kinematic (Lorentz) 7 
boost, which is especially important in the out direction. 
The event generator Therminator 



14, LL5fl can shed 



more light on the source breakup and emission dynamics. 
It gives thermal emission from a longitudinally oriented 
cylinder of radius /O max , includes all known resonance de- 
cays, assumes Bjorken longitudinal boost invariance and 
Blast- Wave transverse expansion with radial velocity v r 
semi-linear in p [lfj|, i.e. v r (p) = (p/p ma , x )/(p/p max + v t ), 
where v t = 1.41. A differential fluid element is a ring 
defined by cylindrical coordinates z and p\ it breaks up 
at proper time r in its rest frame or at time t in the lab 
frame, where t 2 = t 2 + z 2 . The freeze-out hypersurface 
is given by r = To + ap, where tq is the proper breakup 
time for p =0 and a is the space-time correlation parame- 
ter. In Blast- Wave mode, Therminator sets a ~ —0.5 for 
source emission or "burning" from outside in as in many 
hydrodynamical models. 

Using a set of parameters tuned to fit charged pion 
and kaon spectra 17| , pion pairs from Therminator were 
obtained with the effects of all known resonance decay 
processes on and off. These pairs were then transformed 
to the PCMS, following the same set of kinematic cuts 
and Lorentz transformations as in the data analysis, to 
obtain S(r^) distributions for comparison with the data. 

Figure [3] shows that the 3D source function gener- 



ated by Therminator calculations (solid triangles) with 
To = 8.55 fm/c, p m ax = 8.92 fm and other previously 
tuned parameters (l7j , underestimates S(r x ), S(r y ) and 
S(r z ). Open triangles (Fig. [3]) show that resonance de- 
cays reproduce S(r y ) ((b)) and extend the calculated 
source function in x ((a)) as expected, but not enough to 
account for the long tails in S(r x ) and S(r z ). This sug- 
gests that the latter have substantial contribution from 
pion pairs with significantly longer emission time differ- 
ences. Attempts to fit the distributions by only increas- 
ing To or with a > failed, suggesting a fireball burning 
from outside in. 

The generated distribution of time differences can also 
be lengthened by sampling pions from a family of hyper- 
surfaces defined by a range of values of proper breakup 
times t'. One such parametrization consists of replacing 
t by t' chosen from an exponential distribution dN / dr' — 
e( -^ T r ^ exp[— (t' — t)/At]. In this parametrization, the 
width of the distribution At represents the mean proper 
emission duration in the rest frame. Figure [3] shows that 
this approach, with At = 2 fm/c (open circles), leads to 
a fairly good match to the observed source profiles in all 
three directions. A 10% change in At spoils this match. 

Figure [3J^d) shows the relative emission time distribu- 
tion in the LCMS, AtLCM, for pion pairs from events 
with the parameterizations indicated. For a fixed To = 
8.55 fm/c (At = 0) and resonance decays excluded, the 
distribution AtLCM is narrow, (|AtLCM|) = 2.4 fm/c. 
The addition of resonance decays adds a long tail and 
gives (|AtLCivi|) = 8.8 fm/c. Replacing t with the ex- 
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ponential distribution r' with At = 2 fm/c, results in a 
AtLCM distribution which is significantly broadened to 
give (|AtLCM|) = H-8 fm/c. The wider distribution of 
time delays is needed to reproduce the source distribu- 
tions. This implies a finite non-zero proper emission du- 
ration in the emission rest frame. Note that this AtLCM 
distribution broadening has only a small effect on S(r y ). 

Figure[3]shows that substantial time differences AtLCM 
are required by the source distensions; however, the in- 
terplay between proper time and breakup dynamics is 
model dependent. The picture which emerges from the 
data, in the context of the Therminator model, is con- 
sistent with that of an expanding fireball with proper 
breakup time To ~ 9 fm/c, which hadronizes and emits 
particles over a short but non-zero mean proper emission 
duration At = 2 fm/c. 

In summary, a new model-independent, three- 
dimensional source imaging technique has been applied 
to extract the 3D pion emission source function in 
the PCMS frame from Au+Au collisions at ^/snn = 
200 GeV. The source function has a much greater extent 
in the out (x) and long (z), than in the side (y) direction. 
Therminator model comparison indicates a fireball burn- 
ing from outside in with proper lifetime To ~ 9 fm/c and 
a mean emission duration At ~ 2 fm/c, leading to sig- 
nificant relative emission times ((|AtLCM|) ~ 12 fm/c), 
including those due to resonance decay. 
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